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Abstract 

Extending the Standard Model (SM) with sterile (“right-handed”) neutrinos is one of 
the best motivated ways to account for the observed neutrino masses. We discuss the 
expected sensitivity of future lepton collider experiments for probing such extensions. 
An interesting testable scenario is given by “symmetry protected seesaw models”, which 
theoretically allow for sterile neutrino masses around the electroweak scale with up to 
order one mixings with the active (SM) neutrinos. In addition to indirect tests, e.g. via 
electroweak precision observables, sterile neutrinos with masses around the electroweak 
scale can also be probed by direct searches, e.g. via sterile neutrino decays at the Z 
pole, deviations from the SM cross section for four lepton final states at and beyond 
the WW threshold and via Higgs boson decays. We study the present bounds on sterile 
neutrino properties from LEP and LHC as well as the expected sensitivities of possible 
future lepton colliders such as ILC, CEPC and FCC-ee (TLEP). 
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1 Introduction 


The origin of the observed neutrino masses is one of the great open questions in particle 
physics. There are various ways to introduce massive neutrinos, which all require an ex¬ 
tension of the particle content of the Standard Model (SM), or the introduction of effective 
operators which have to be generated at some higher energy scale by new physics involv¬ 
ing additional particles. Currently, no experimental evidence exists to select between the 
various proposed extensions of the SM towards massive neutrinos. 

One of the best-motivated and most minimal extensions of the SM for providing neu¬ 
trino masses consists in adding “right-handed” (often named “sterile”) neutrinos to the 
SM degrees of freedom. Among the types of fermions within the SM, i.e. up-type quarks, 
down-type quarks, charged leptons and neutrinos, the neutrinos are the only type without 
a right-chiral counterpart. In Grand Unified Theories (GUTs) based on the gauge group 
SO(IO), for instance, which has a left-right symmetric particle content, the “right-handed” 
neutrinos are therefore predicted. 

“Right-handed” neutrinos would be SM-gauge singlets, and because of this they are 
often referred to as “sterile”. Nevertheless, they can interact with the SM particles via 
their Yukawa couplings to the lepton doublets and the Higgs doublet. This coupling results 
in a Dirac-type mass for the neutrinos when the neutral component of the Higgs develops 
a non-zero vacuum expectation value. Furthermore, as gauge singlets, the sterile neutrinos 
can also have a mass term with their charge conjugates, i.e. a Majorana mass term. This 
leads to a mixing between the active and the sterile neutrinos. In the mass basis both, the 
light and the heavy eigenstates, couple to the Z and the W bosons. 

With n sterile neutrinos, the full neutrino mass matrix would be a (3 -|- n) x (3 -|- n) 
matrix. Out of the 3 -|- n mass eigenstates, at least three have to be light, i.e. below about 
0.5 eV, in order to account for neutrinos oscillations, cosmological observations as well as 
constraints from neutrinoless double beta decay. The masses of the other (mainly sterile) 
neutrinos, as well as their couplings to the SM particles, are basically free parameters of 
the theory. 

Sterile neutrinos in various mass ranges have been discussed in the literature (see e.g. [1] 
for a review): For instance, sterile neutrinos with masses in the eV range could lead to 
effects in short distance neutrino oscillation experiments by introducing an additional mass 
squared difference. keV mass sterile neutrinos are candidates for “warm” dark matter, and 
very heavy sterile neutrinos around Mgut ~ 10^® GeV are often predicted from GUTs. 

Here, we are interested in sterile neutrinos with masses around the electroweak (EW) 
scale, such that direct searches at present and future colliders are possible. Gompared to our 
recent work [2] where we assumed the sterile neutrinos to have masses sufficiently above the 
EW scale to test them via probes of non-unitarity of the effective leptonic mixing matrix, 
sterile neutrinos with EW masses can now be produced on-shell in particle collisions. In 
addition to indirect tests, e.g. via electroweak precision observables, they can now also be 
tested via, e.g., sterile neutrino decays at the Z pole, deviations from the SM cross section 
for four lepton final states at and beyond the W pair production threshold, and via Higgs 
boson decays. 

In this paper we first study the present constraints on sterile neutrino properties in this 
mass range, including the whole relevant data from indirect tests (as in m as well as the 
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present constraints from LEP and LHC on the above mentioned processes. Furthermore, 
we will provide first estimates for the expected sensitivities of future colliders such as the 
International Linear Collider (ILC), the Circular Electron Positron Collider (CEPC) and 
the electron-positron mode of the Future Circular Collider (FCC-ee, formerly known as 
TLEP), for testing EW sterile neutrinos and compare the prospects of direct and indirect 
searches. 

The paper is organized as follows: Section 2 contains a description of a minimal symme¬ 
try protected “low scale” type-I seesaw scenario, which we will use as benchmark framework 
for our analysis. In section 3 we derive the constraints from the present data, including 
“direct” and “indirect” searches for sterile neutrinos. In section 4 we present estimates for 
the sensitivity of future lepton colliders to sterile neutrino properties. Section 5 contains a 
discussion of our results and the conclusions. 

2 A symmetry protected “low scale” type-I seesaw scenario 

As described in the introduction, we investigate sterile neutrinos with masses around the 
EW scale. Such sterile neutrino masses can be realized in a “natural” way together with 
large (even 0(1)) Yukawa couplings to the lepton doublets and the Higgs doublet if there is 
a “lepton-number-like” symmetry which controls the size of the light neutrinos’ masses, i.e. 
protects them from getting too large 0 In fact neutrino masses in such scenarios are small 
when this protective symmetry is only slightly broken, in contrast to the usual seesaw 
mechanism where the smallness of the light neutrinos’ masses comes from the heaviness 
of the sterile states. In our analysis we will focus on a minimal version of the symmetry 
protected scenario, where the experimentally observable effects stem from one pair of sterile 
neutrinos (having opposite charges under the protective symmetry). 

2.1 The scenario: Extension of the SM by EW scale sterile neutrinos 

To realize a low scale seesaw mechanism with a pair of sterile neutrinos (/ = 1,2) 
without highly suppressed neutrino Yukawa couplings, we impose e.g. a “lepton-number- 
like” (global) U{1) symmetry, where has the same (opposite) charge as the left- 

handed SU{2)l doublets L", a = e,^,T. Neutrino masses arise when this symmetry gets 
slightly broken, as e.g. in the so-called “inverse” mM or “linear” [5] variants of the type I 
seesaw mechanism (see also e.g. [SHS]). The Lagrangian density in the symmetric limit is 
given by 

^ = ^SM-^M+ B.c. , (1) 

where contains the usual SM field content and with L“ and cj) being the lepton and 
Higgs doublets, respectively^ The y,y^ are the complex-valued neutrino Yukawa couplings 
and the sterile neutrino mass parameter M can be chosen real without loss of generality. 

^The term “natural” is understood here in the’t Hooft sense: Setting the masses of the light neutrinos 
to zero enhances the symmetry of the theory. 

■^We remark that the “lepton-number-like” symmetry mentioned above is just an example and basically 
any symmetry leading to the above Lagrangian (with effects from possible additional terms being sufficiently 
suppressed) may be used. 
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A third (or even more) sterile neutrinos may exist in addition, but we assume that 
it has (or they have) zero charge under the “lepton-number-like” symmetry such that in 
the symmetry limit they decouple from the other particles (since no Yukawa couplings to 
the lepton doublets are allowed and they also cannot mix with the other sterile states) 
and will be ignored. In this case, are the relevant parameters for studying 

the phenomenological consequences of a pair of EW scale sterile neutrinos H Only when 
the protective symmetry gets broken, all sterile neutrinos contribute to the generation 
of the light neutrino masses and all three light neutrinos will obtain small masses. For 
phenomenological tests of the scenario at colliders or low energy precision experiments, 
the very small symmetry-breaking terms have negligible effects (see e.g. mm and we will 
therefore study the limit where the symmetry is intact. 

After electroweak symmetry breaking, one obtains the following 5x5 mass matrix of 
the (relevant) electrically neutral leptons: 
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+ H.C., 


( 2 ) 


with the Dirac masses ma = Vva'^'ESN !where y^^ are complex coupling constants and 
uew = 246.22 GeV. Note that in this limit of exact symmetry the right-handed neutrino 
does not couple to the SM leptons, and that the three lightest neutrinos are forced to 
be exactly massless. Diagonalising the mass matrix from eq. ([2]), which we will denote by 
A4, with the unitary matrix U yields the mass eigenstates: 


MU = Diag (0,0, 0, M, M) , 


(3) 


where we have neglected 0{6‘^) corrections to the masses of the heavy neutrinos. The 
(complex) active-sterile mixing parameters are defined as 


_ ytg VEw 
~ V2 M 


(4) 


and the quantity 

= ( 5 ) 

a 

such that to second order in the mixing parameter 6^ the leptonic mixing matrix U is 

® We note that in the specific case that indeed no further sterile neutrinos exist, i.e. only the two which 
form the pseudo-Dirac pair, then the lightest neutrino remains massless and there are correlations between 
the , y^,^ which depend on the elements of the light neutrino mixing matrix [S]. Similarly, such 

correlations also arise in a special limit of the R-matrix parmeterisation [9]. In this study, we will not 
impose additional constraints of this type but consider the , Vv^ as independent parameters. 


3 











yua 

9a 


Vva ~ 

- 

y/2M 


9q = 

f’EW 

- 


^aj3 — 

_ 2M^ 

-0*a0y 

- 


Table 1: Relation between the sterile neutrino parameters and 9a, and the leptonic non-unitarity 
parameters Sap used in ref. [2]. 
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( 6 ) 


The elements of the non-unitary Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix M 
can be written as [TO] 

{u^)pi , (7) 

with Ui being a unitary 3x3 matrix. 


2.2 Parameters 


As described above, the relevant new parameters of the scenario under consideration are the 
three complex Yukawa couplings yu^,yufj,,yuT and the mass M. Via eq. dH), these parameters 
can be mapped onto the three (also complex) active-sterile mixing parameters 9e,0^,9r- 
Concerning physical processes, where the sterile neutrinos are very heavy compared to 
the experimental energy scale, they can be integrated out. The remaining effect is then 
given by the lepton-number conserving dimension six operator with coefficients 


d=& _ 

m2 ’ 


( 8 ) 


which causes an effective non-unitarity of the leptonic mixing matrix via a contribution to 
the kinetic terms of the active neutrinos, as discussed e.g. in [2|. The deviation of the PMNS 
matrix from unitarity, i.e. Sap = is obtained either from the coefficients in 

eq. (|8|) or from the definition of the mixing matrix J\f in eqs. (|6|) and d?]) directly. To leading 
order in the mixing parameters the relation is given by 




?;2 

^EW 

2 


= -9*j0 , 


(9) 


with the definition of the mixing 9a from eq. O- We summarise the parameters and the 
relevant mappings in tab. [TJ 
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2.3 Modification of the weak currents 


We can collect the left-handed neutrinos and the charge conjugate right-handed fields 
{N^y into the column 

n= ■ ( 10 ) 

The mass eigenstates are given as 

T 1" 

flj = {vi, U2, 1 ^ 3 , ^4, N5)j = Uj^Ua ■ ( 11 ) 

Now we can write down the weak currents in the mass basis: 

5 

Pl Uai fli + H.c. , (12) 

i=l a=e,fi,T ' 

5 

jfj. = nj 7^ Pl Uai fii , (13) 

. . z cw 

2 j = l Q;=e,/i,r 

where g is the weak coupling constant, cw is the cosine of the Weinberg angle and Pl = 
^(1 — 7 ^) is the left-chiral projection operator. It is convenient to dehne the quantity 

E (14) 

a=e,li,T 

such that with i < 3 and j = 4, 5 we cobtain 

0'=e,fi^T ^ (y.=e^ii^T ^ 

where the non-unitary PMNS matrix M was defined in eq. d?]). The gauge couplings of 
the fermion current with two heavy neutrinos are proportional to 'djk for j, k = 4,5, which 
satisfy 

\«,t\=y\ (16) 

with 0^ defined in eq. ([5|). With the above definitions, the weak currents involving the 
heavy neutrinos in the mass basis can be written compactly as 

D ^eaia 7 ^.PL{-m + Ny + }^.c., ( 17 ) 

5 

E . (18) 

^ C\Y . . 

^.7=1 


Analogously we can express the Yukawa part of the Lagrangian density in the mass basis, 

E + H.c. ^Ey'^’-E / Pc^ini + H.c. 




V2M 

ve-w 


j=4,5 


2=1 


+ H.C., 


(19) 


with being cjy the neutral component of the Higgs doublet. 
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2.4 Decay rates involving heavy sterile neutrinos 

With the weak currents in eq. (ED and (fTBp and the Yukawa terms in eq. (fTni) . the heavy 
neutrinos N 4 and Y 5 couple to the weak gauge bosons and the Higgs boson, respectively. 
They can either be produced in decays from gauge and Higgs bosons, or decay into leptons 
and bosons, depending on which process is kinematically allowed. 

First, we consider the case of M < mu/j "iz, "T/d which yields the following decay 
channels: 


W+^£N, W-^IN, Z^Nv, Z^NN, h^Nv. ( 20 ) 

We have suppressed here the indices of the neutrino mass eigenstates and of the leptons. 
VF, Z are the weak gauge bosons, h is the SM Higgs boson and £ = e, fj,, t denote the 
charged leptons. Note, that for the Z and Higgs boson decays also the Hermitean conjugate 
processes have to be taken into account. Neglecting the masses of the light neutrinos and 
charged leptons the corresponding decay rates for i = 1 , 2 ,3, j, A: = 4, 5 are given by 

| 6 » q ,|2 Gf ml 


r(w - 
r(w+ - 

T{Z 
T{Z - 
T{h 


Njia, 


Njia 


‘-W 


VI 


NjNk) = \^jk\ 


2 Gy/ 2 TT 
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11 ( 1 + 1 ) 


I#- '2 


Gp m^ 


vr 


6^2' 
,Gf m^ 


11(1+1) (/^Z) : 




6 V 2 Tr 
mh I'&ij 






^11(2)(a^z) , 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


where we introduced fix = M/mx, Gp the Fermi constant, and the kinematic factors 


are 


11(i+i)(a^x) — 2 “ 1^^) (2 + /^^)’ 

n(2)(iijv) = 2 ~ + 


(26) 

(27) 


We note that the decay rates of Z to Y are the same as the ones to N, e.g. F(Z —)• ViNj) = 
T{Z —)• ViNj). Both processes have to be taken into account when calculating Rinv^ as will 
be discussed below. 

To obtain the total Higgs decay rate into neutrinos, we observe that to leading order in 
the active-sterile mixing parameters 

3 

= l^el^ + ■ (28) 

i=l 

Therefore, by summing eq. (I25p over j = 4, 5 and i = 1,2,3, and including the Hermitian 
conjugate process, we obtain to leading order in the mixing parameters 

mh 0 ^ 


T{h vN) = 


2\2 




(1 - Ri) 


(29) 
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The complementary processes which are kinematically available for M > 


namely the corresponding decay rates for the heavy neutrinos, i.e. j = 4, 5, are 

Nj W ia, ^ ■ ( 30 ) 

The corresponding decay rates for i = 1,2, 3, j = 4,5 are given by 

r{N,^W£a) = ( 31 ) 

r{Nj^Zn) = (32) 

TiNj^hn) = (33) 

We use the following parametric values m- 


Parameter 

mz [GeV] my/ [GeV] 

rrih [GeV] 

Gp [GeV-^] 

Value 

91.1875 

80.358 

126.0 

1.1663787x10-^ 


3 Present Constraints 

Before we study the sensitivities of future colliders, we discuss the constraints on sterile 
neutrino properties from the currently available experimental data. We start with “indirect” 
constraints from precision tests of the SM and then turn to “direct” tests focusing on sterile 
neutrino decays at the Z pole, deviations from the SM cross section for four lepton final 
states at and beyond the WW threshold and Higgs boson decays. 

3.1 “Indirect” constraints from precision tests of the SM 

We consider the mass of the heavy neutrinos, M (note that we have only one mass scale 
here due to the protective symmetry), to be in the range from ~ 10 GeV to ~ 250 GeV. In 
the presence of the sterile neutrinos the theory predictions for various precision observables 
get modified. In this subsection we extend the analysis of our recent work [2], where we 
assumed the heavy neutrinos to have masses sufficiently above the EW scale to test them 
via probes of non-unitarity of the effective leptonic mixing matrix, to masses around the 
EW scale. A discussion of constraints on sterile neutrinos in the mass range below 10 GeV 
can be found in ref. m- 

Eor various observables, where the experiments are performed at energies much below 
M, the results from [2] still apply (since the effective theory treatment is still applicable) 
and we can simply translate them into constraints on the sterile mixing parameters using 
table [TJ We will mainly revisit the observables where there are changes due to M around 
the EW scale, such as the electroweak precision observables (EWPOs) at colliders as well 
as the rare charged lepton flavour violating (LEV) decays. 
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Effects on the Fermi constant 


The Fermi constant Gp is measured from muon decays which get modified due to the effects 
of the sterile neutrinos in the charged current interactions. Denoting the Fermi constant 
extracted from muon decays as G^, we obtain the relation (at tree-level and to leading 
order in the mixing parameters) 

Gl = GUi - |0ep)(l - \0^\^) ■ (34) 

This has consequences for the theory predictions of many precision observables. 


Electroweak precision observables 

The tree level relation between sw = sin^Wj Gp and a is given by 

2 2 _ a(mz)vr 

— 


V2Gpm?^ ’ 


(35) 


which yields the following theory prediction for which is modified with respect to the 
SM due to eq. (l34l) : 


2 


1- Wl- 


2 ^ 2 , 


an 






(36) 


Furthermore, together with the tree-level relation m\(^ = the theory prediction for 
the W boson mass is modified to 


(37) 


I 

with the weak mixing angle syi/ from eq. (I36p . 

Z boson decay parameters 

The modification of the Fermi constant in eq. (I34|) also modifies the tree level decay rate of 
the Z boson into fermions //. For f ^ p we have 




(38) 


6V27r ^(1-|0^|2)(1-|0,12)’ 
with Nc being the colour factor and gvj-, dAj the vector and axial vector coupling 

9VJ = tI - 2QfSw , 9A,f = tI , (39) 


with the third component of the isospin and the electric charge Qf. Through G^ and 
Sw-, the decay rate Fjj is affected by the modification of the light neutrino couplings. An 
observable which is very sensitive to the modifications due to the mixing parameters is the 














decay rate of the Z boson into two light neutrinos, i.e. with i, j <3 and to leading order in 
the active-sterile mixing parameters: 

3 

- Khf X rz^-,SM X [(1 - 0e|')(l - \0^\^)] " , (40) 

i,j=l O',13 


where = G^m|/(6\/27r) is the decay width for Z —)• Uai^a for a specific flavour, in 

the SM. 

We now turn to the hadronic pole cross section and the invisible decay rate Rinv 
of the Z boson, defined as 


^0 _ 127r Yz^eeF Z^had jy _ ^Z—^inv 

^had — XfJ U2 ’ ~ “r 

^-‘z ^ Z ^ 


(41) 


where Tz^inv is the invisible partial decay width and Tz^had is the sum over all the 
hadronic partial decay widths of the Z boson and is the total Z decay width. To 
leading order in 9q. we obtain the following expression 


^ - 0.09 

where we used the values for the parameters = —0.82 and cr = —0.67. If the heavy 
neutrinos would not decay inside the detector, the would be —1. For estimating the 
parameters Co-,r we have assumed that the heavy neutrinos decay within the detector (as will 
be the case for most of the considered parameter space) and that the kinematically available 
SM fermions are massless. Furthermore, following m, we assumed that all processes where 
the heavy neutrinos decay into hadrons (and a light neutrino) are counted as hadronic 
events, whereas all leptonic and semileptonic N decays are rejected by the event selection 
hlters. This estimate is sufficient for the discussion in this paper, but should be replaced 
by a more accurate treatment at latest when a signal is found. 

It is often stated that for heavy neutrino masses much smaller than niz unitarity of 
the PMNS matrix is effectively recovered and that therefore the prediction for the invisible 
Z decay rate coincides again with the SM one (with = 3). However, as long as the 
heavy neutrino is heavier than the muon, there is in any case a dependency on the mixing 
parameter combination + due to the use of the Fermi constant as input parameter. 

Analogously to the above discussions for Rinv and we also include the pseudo¬ 

observables Ri,Rb and Rc- 

Ri = [i?dsM(l + O.15(l^e|' + |0Ml")-O.O702n(i+i)(/iz)), (44) 

Rb = [i?fe]sM(l-O.O3(|0e|" + |0Ml")-O-OOl^'n(i+i)(/^^))> (45) 

Rc = [i?c]sM(l + O.O6(|0e|" + |0/.P)-O.OOO302n(i+i)(/iz)). (46) 

We assume that the sterile neutrino decays do not signihcantly affect the experimental de¬ 
termination of niw and the effective weak mixing angle (such that they are only sensitive to 
+ due to Gi? as given in eqs. p6]l and ([37)1 ). The SM predictions and experimental 
values for the EWPOs are taken from ref. |14] . 


(|0eP + |0MP) > (43) 


Ri 


[-Rm,;]sM + CH 71(1+1) if^z)) 
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r)E 

Ku 

T)'K 

DTT 

Ku 

R^e 

Process 

r(r^i>TePe) 

r(fl^Up,ePe) 

r(TT—>-ePe) 

ri-K^iiPu) 

r{T^KVr) 
r{/f— vfiDu) 

T(T^KVr) 

r(K^ePe) 

Bound 

1.0018(14) 

1.0006(21) 

1.0021(16) 

0.9956(31) 

0.9852(72) 

1.018(42) 


Table 2: Tests of lepton universality used in our global fit. Experimental results are taken from ref. m- 


Lepton universality observables 

The lepton universality observables considered here are defined as ratios of decay rates: 

= Tq /r^, where denotes a decay width including a charged lepton Iq and a neu¬ 
trino. They are defined such that in the SM = 1 holds for all a, j3, X, as a consequence 
of lepton universality. We include here constraints from vr, /t, r and K decays which stem 
from experiments at comparatively low energy and which are currently dominating the 
constraints. This allows us to use the results from [2] for these processes, translating the 
parameters using table [TJ Constraints on sterile neutrinos from lepton universality tests 
have also been studied recently in refs. |15lll6j . We note that in contrast to [2] we are not 
including W decays here, which however only had a negligible impact on results of the fit 
in [2]. The active-sterile mixing between the left-handed and right-handed neutrinos leads 
to modified theory predictions of the form 

= 1 - ^ (l««l" - \h\^) . (v) 

which thus allows to probe differences between the Oa- We display the present experimental 
constraints on the universality observables in tab. [2j 


Rare flavour-violating charged lepton decays 

The decay rate for lepton flavour violating charge lepton decays Ip £o-7 are given by 




20487r4 




k=l 


(48) 


where terms ~ are neglected and where F{xk) is a loop-function which 

depends on the mass ratio Xk = \mu^/M\Y\\ 


F{x) 


10 — 43x -b 78x^ — 49x^ -b 4x^ -b 18x^ In x 
3(1 -x)4 



(49) 


The nii/f, are the mass eigenvalues of the light and heavy neutrinos, and we shall approximate 
them with 0 and M, respectively. By using the unitarity of the neutrino mixing matrix U 
(up to second order in 9a), we can write 


E V,kUlF(xt) 


k=l 


3 5 

E V,tUlF(0) + E U,uUIf(xm) 


k=l 


k=A 


-ep9l[F{Q)-F{xM)] , 


(50) 
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Process 

Prediction with active-sterile neutrino mixing 

90 % C.L. bound 

Br{y eq) 
Br{T —>• eq) 
Br{T —>■ /iy) 

2.4 X 10-^ (1 - 0.3F(xm)) \9p9^\'^ 

4.3 X 10-4 (1 - 0.3 F{xm)) \0rde? 

4.1 X 1O-4(1-O.3F(xm))|0*0„P 

5.7 xlO-4^ 

1.5 xlO"® 

1.8 xlO"® 


Table 3: Present bounds on the charged lepton flavour violating processes la —>■ ^/37 and predictions in the 
presence of sterile neutrinos. The experimental bounds on /r —>■ 67 are from the MEG collaboration |18| . 
the ones on r decays are taken from ref. m- The function F{x) is defined in eq. (1491) . 


with xm = M/mw- Plugging this into eq. (H8|) . the branching ratio for the process £p —>■ 
can be expressed as 


Brp„ = Br{ip 


_ 100a 2i/i 

^'pZ^fT^cr) “XT \Sp\ \0o 
OOTT 


1 - y^F{xm) 


(51) 


Notice, that with M —0 unitarity is restored, as it should be. Furthermore, we remark 
that in the limit M S> the above expression seems to differ from the corresponding 
one in ref. [2], where the low energy effective theory was considered. Note however that 
“decoupling” automatically implies 0 —)• 0, such that in both frameworks the low energy 
effects of the extra sterile neutrinos disappear as they should. The present experimental 
constraints on rare charged lepton decays are summarised in tab. [3l 


Other precision constraints 

In addition to the “indirect” tests mentioned above, we also include the constraints from 
the NuTeV experiment, from CKM unitarity tests and from low energy measurements of 
s^. To calculate these additional constraints we follow the analysis of [2], translating 
the parameters using table [TJ The modified theory predictions for the EWPOs, NuTeV 
observables and low energy measurements of the weak mixing angle due to active-sterile 
neutrino mixing are summarized in tab. 01 The constraints from CKM unitarity tests are 
discussed in ref. [2] in sec. 3.2.3. 

3.1.1 Present constraints from “indirect” tests: Global fit resnlts 

To obtain the present constraints on sterile neutrino extensions of the SM from precision 
observables, we perform a Markov Chain Monte Carlo (MCMC) fit for the three Yukawa 
couplings and mixing parameters, and 9a-, for 10 GeV < M < 250 GeV, and extract the 
highest posterior probability density (HPD) intervals at 90% confidence level (CL). We use 
the experimental constraints discussed in the previous section, which are essentially based 
on the observables also used in ref. [2] unless stated otherwise in the text, adapted to sterile 
neutrino scenarios with M in the EW range. 

We display the resulting upper bounds for the parameters \9a\ and \yua\^ with a = e, y, r, 
in fig. [H For comparison, we also show the exclusion limit from the direct searches of the 
LEP-I experiment Delphi, which will be discussed below. 

We find, in agreement with the results in refs. [2l[23], that for {yuj there is also a non¬ 
zero lower bound at 90% confidence level. In the following, however, we will only use the 
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Prediction with heavy neutrinos 

Prediction in the SM 

Experiment 

[wiv]sm(1 + 1'^ + 

[r'leptlsM{l + 0-59{|(?ef + ^}) 

[{»'A(r)"lsM{l - 0.71(19.1" + |9„|")) 
[{4 “)"Ism{1 - 0.71(19.1" + |9„|")) 

80.358(8) GeV 
83.966(12) MeV 
0.23150(1) 
0.23150(1) 

80.385(15) GeV 
83.984(86) MeV 
0.23113(21) 
0.23222(27) 


ii,]sM(l + O.3|0e|"-1.7|0^|^) 
i?p]sM(l + O.l|0ep-1.9|0^|2) 

0.3950(3) 

0.4066(4) 

0.3933(15) 

0.4034(28) 


+ 9-l( 

(1 + 15.1( 0eP + 

-73.20(35) 

0.0710(7) 
1.520(24) X10"'^ 

-72.06(44) 

0.064(12) 

1.31(17)xl0-^ 


Table 4: Experimental results, SM predictions and the modification in the presence of sterile neutrinos for 
mw, the effective weak mixing angle, the NuTeV observables and for the low energy measurements of the 
weak mixing angle. The SM predictions and experimental values for the EWPOs are taken from ref. |14| . 
The values of taken from Ref. [^. The NuTeV results on deep inelastic scattering 

of neutrinos and anti-neutrinos on nuclear matter has been taken from ref. m- The theory uncertainty 
stems from s^. The results on are from Ref. [2^. For we used Sw{Mz) ~ 0.2315, and its 

error is dominated by the uncertainty of the radiative QED correction factors. 


upper bound on as constraint. The best fit value for is zero and the uncertainty 
on \yu.r\ is much larger than that of the other two parameters. 

3.2 Present constraints from “direct” searches 

We now turn to the current constraints from “direct” searches, i.e. via sterile neutrino decays 
at the Z pole, deviations from the SM cross section for four lepton final states at and beyond 
the WW threshold, and via Higgs boson decays. As we will discuss below, especially the 
first two tests provide the strongest constraints for specific mass ranges (below M ~ 150 
GeV), whereas the “indirect” tests are more sensitive for larger M. In the next section 
we will estimate the sensitivity improvements which could be possible at envisioned future 
colliders. 

3.2.1 Search for sterile neutrinos produced in Z boson decays 

The LEP-I collaborations Delphi [23], Opal [25], Aleph [26] and L3 [27] have performed 
analyses searching for “heavy neutral leptons” - or sterile neutrinos - at the Z-pole. The 
Feynman-diagram for sterile neutrino production which is dominant at the Z pole is shown 
in fig. [21 The results of the LEP collaborations can be expressed as an upper limit on the 
branching ratio for Z bosons decaying into a light and a heavy neutrino. It can be used to 
constrain the sterile neutrino parameters as we now discuss: 

The strongest bound on the branching ratio for the processes Z ^ v N comes from the 
Delphi collaboration. It is given at 95% C.L. as 

^ < 1.3 X 10"®, (52) 

which includes the processes Z i>j <3 Nj, j = 4,5 and the Hermitean conjugate processes. 
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Figure 1: Upper limits on the active-sterile neutrino mixing parameters from “indirect” tests at 90% CL. 
The left panel shows the mixing parameters Oa, a = e,fj,, r, the right panel shows the Yukawa couplings . 
The purple line represents the direct search constraints on the parameter space from Delphi |24| . 







Figure 2: Feynman diagram dominating the production of sterile neutrinos at the Z pole. 
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Figure 3: Feynman diagram dominating the production of sterile neutrinos at the WW threshold. 


With the expression for the corresponding decay rate in eq. (1231) . the experimental upper 
bound from eq. (I52h can be used to put upper bounds on the sum over all the active-sterile 
mixing parameters: 


< 1-1 X 10 


-5 


(l-/r2)2(2 + ;,2) 


(53) 


The resulting constraint is shown in figUl 


3.2.2 Search for sterile neutrinos in 4£ final states at LEP-II 

At LEP-II, the properties of W bosons were studied at center of mass energies -y/i at and 
beyond the WW threshold. One of the relevant observables is the cross section for WW 
production, which can be reconstructed from the W decays into four-lepton hnal states, 

e+e“ ^ WW P£-£+i^ . (54) 

The observed cross section was found to agree with the SM prediction and the Aleph 
experiment at LEP-II has placed a bound on possible SM deviations dAieph, dehned via 
|(5cJu/rF—>.4£| ^ <5exp at IfJ C.L. [28]: 

Aleph 

sAleph = 1 - = 0.005 ± 0.011,tai ± 0.007,y,t. (55) 

™VF1F 

For our analysis we will combine the two contributing uncertainties in quadrature. For the 
SM prediction, we use the cross section for WW production from RacoonWW [29] . 

The dominant correction caused by sterile neutrinos to e^e“ —)• 4£ (at the considered 
energies) arises from diagrams of the type shown in fig. (3] which produce the same final 
states from N decays 

e+e- ^vN ^ v£-W^ v£-£+v , (56) 

with a cross section cri,pf^ 4 £. The produced four lepton final states would have been misin¬ 
terpreted as as a contribution to WW pair production, which allows to constrain the sterile 
neutrino properties from the Aleph bound of eq. (I55p . In the narrow width approximation, 
we obtain: 

i'max 

<TuN^ai= J (57) 

^min 
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Channel 


Rww 

Rzz 

Atlas 

-1 1 y+0.27 

-‘^■-‘^'-0.27 

1.08t(!;^2o 

1 44+0-40 
4-^^-0.33 

CMS 

1 14+<J.3U 

n 70 +U. 2 U 

U./Z Q 18 

U.yo 0.25 

combined 

1.15(27) 

0 . 88 ( 20 ) 

1.11(30) 


Table 5: Currently best measured decay ratios Rxx = Br{h —>■ /Br(h —> from CMS |32H34| 

and ATLAS [35]. 


with Br[N —>■ 3£) ~ 0.2 (not counting decays to neutrinos) and the kinematic limits 
train = ~'S/2(1 + /3) + 771^ and tmax = ~'S/2(1 — /3) + m^. The differential cross section is 
given by [50] 


da, 


e+e~ ^uN 

dt 


/-i2 4 '5 

2 TT 5^ 

2=1 


2|t9j4|^t^ 


(s — m|)2 


{t{t — M^) + u(u — M^)) 


+ 


+ 


(1 ~ 2 S^) 




(1 ~ 2 


t-m^ 


2(s — m^) u — m 


w 


u(u — M^) 
t (t - M^) 


(58) 


where swi cw,t\Y are the sine, cosine and tangens of the weak mixing angle 9w, respectively. 
The bound on \9e\ is then obtained from the requirement 

— ^exp 5 (^9) 


which is shown below in fig. [9| by the black line and the grey area and in the summary plot 

fig. [m 

A different analysis searching for sterile neutrinos beyond the Z mass threshold has 
been conducted by the L3 collaboration m- They consider the production of active and 
sterile neutrino and the subsequent decay chain N igW iej{j), where j is a hadronic 
jet. The reconstructed invariant mass of the heavy neutrino would manifest as a peak in the 
invariant mass distribution. The bounds of |3T] are of the same order as the ones presented 
here. 


3.2.3 Higgs boson decays at the LHC 


We now consider the constraints from the present LHC measurements of the Higgs decay 
parameters, which are shown in tab. |5l Due to the large uncertainties on the fermionic 
branching ratios, we focus here on the decays h —>■ VV, for the vector bosons V = 7 , Z, W. 
For M < mh, the heavy neutrinos are produced in Higgs decays and will modify the 
branching ratios, which allows to constrain the neutrino Yukawa couplings (or equivalently 
the heavy-light mixing angles). There are basically two effects of the additional Higgs 
decays: 

Firstly, the total Higgs decay width is enlarged, which effectively reduces all the SM 
branching ratios by a factor 


r/i,SM + ^h^uN 


(60) 
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Figure 4: Upper bound on the branching ratio of the Higgs boson into heavy and light neutrinos from 
“indirect” tests (solid blue line) and the direct search by Delphi at LEP-I [53] (solid purple line). The 
branching ratio Brh^vN ■= j Br{h —>■ ViNj) denotes the sum of the processes h ViNj, for i = 1, 2, 3 
and j = 4, 5 and the Hermitian conjugate processes. For M > mz the present “indirect” constraints allow 
for even an order one Brh^vN- 


In fig. m the upper bound on the branching ratio for the total decay width of the Higgs 
boson into light and a heavy neutrinos is shown by the solid blue line. The plot shows 
that even an 0 ( 1 ) branching ratio would be consistent with the present “indirect” 

constraints. 

Secondly, when heavy neutrino decays take place inside the detector (which is the case 
for the relevant parameter space we are considering here), then the subsequent decays may 
get counted as SM Higgs decays into ZZ or WW since they lead to the same final states. 
Altogether, the experimentally measured branching ratios with respect to the SM prediction 
are given by 


Bvh^xx = r Brh^xx, SM + cxBrh^^x, with cx 


i, X = Z,W 

0, X = 7, / 


( 61 ) 


with Bvh^^x = r/i^i/Ar/(r/i^ sM + Th^-i/Tv)- The experimental precision of the currently best 
measured observables (cf. tab.jSj) can now be translated into sensitivities on the active-sterile 
mixing parameters. 

In fig. El we display the deviation of the Higgs branching ratios into gauge bosons, for 
the squared sum of Yukawa couplings \yuj\‘^ = 10““^, for illustration. Figure [ 6 ] shows 
the sensitivity of the present LHC measurements on the active-sterile mixing parameters 
for the different decay channels, the most sensitive of which is Br{h —)• 77 ). Note that 
when estimating the constraints we assumed that every single heavy neutrino decay is 
counted into the associated decay channel of Higgs to a gauge boson. In reality, some 
events might be removed by the experimental filters. The assumption we made leads to 
a conservative bound, since the two effects described above change the branching ratios 
in opposite directions. Successful filtering would increase the sensitivity of the branching 
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Figure 5: Deviation of the Higgs decay ratios Rxx from the SM prediction Rxx = 1 for the example value 




Figure 6: Constraints on sterile neutrino parameters from Higgs decays at the LHC. 



ratios into W and Z. The estimates obtained here are in good qualitative agreement with 
the analyses in refs. [S] and [3B], where the event signature has been analysed more carefully. 
Other sterile neutrino decay signatures at the LHC have been analysed e.g. in refs. 

[37H1I]. 

4 Possible improvements from future lepton colliders 

We will now estimate how the improved sensitivities of future colliders could allow to test 
sterile neutrino properties. The improvements will concern “indirect” tests, especially via 
the EWPOs, as well as “direct” searches. We will focus on the processes discussed in the 
previous section and consider the ILC, CEPC and FCC-ee (TLEP), which are currently 
discussed, as representative examples. We note already here that the numbers which are 
given are based on estimates for sensitivities as they are currently discussed in the respective 
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Observable 

ILC 

FCC-ee 

CEPC 

CEPC* 

Re 

Rinv 

Rb 

Mw [MeV] 
2,e 

^eff 
al [nb] 

F^ [MeV] 

0.004 

O.OI 

0.0002 

2.5 

1.3 xlO-5 

0.025 

0.042 

0.001 

0.002 

0.00002 

0.5 

1 xlO"® 

0.0025 

0.0042 

0.01 

0.012 

0.00017 

0.5 

2.3 xlO-s 

0.003* 

0.006* 

0.0007* 

0.5 

3.3 xlO"®* 
0.008* 
0.014* 

Reference 

m 

m 

05 ], m 

scaled* 


Table 6: Estimated systematic uncertainties of the ILC, the CEPC and the FCC-ee for future measurements 
of the EWPOs. (The statistical uncertainties would be much smaller.) 

Performance scaled with a factor vTO from FCC-ee, for comparison. See text for details. 


working groups or are extracted from present proposals. These numbers may change and 
affect the comparison between the experiments. 

For instance, ref. [l2] has recently suggested for the CEPC, to increase the integrated 
luminosity off the Z peak, and the use of polarised beams, in order to boost the precision 
of the EWPO measurements. An important result of this suggestion is the improvemed 
precision when measuring syp^e//) which we include in the third column of tabjGl 

The table also includes the current estimates for the possible experimental precision 
supplied by the CEPC study group0. We remark however, that not all information regarding 
the CEPC performance in the electroweak precision sector are available at present. For a 
comparative study of the machine performance we include a second column on the CEPC 
uncertainties in tab. El which contains the precision of the FCC-ee, scaled with a factor 
of vTo. This would correspond to a CEPC with identical performance parameters as the 
FCC-ee but with 10^^ Z bosons instead of 10^^ as considered in [43] for the FCC-ee. 

4.1 Sensitivities of “indirect” searches at future colliders: EWPOs 

As discussed in the previous section, the EWPOs are sensitive to the parameter combina¬ 
tions and |0eP + for a given M. For estimating the possible future sensitvities we 
use the observables and experimental uncertainties listed in tab. El The resulting bounds 
on the parameters \6r\‘^ and l^eP + and the bounds on the Yukawa couplings \yuT-\ and 
^/\yi^e\‘^ + shown in the four plots of £g. [71 Note the relaxed constraint on \6t\ 

(and \yuT-\) for M < mz, due to the phase-space factors in the theoretical prediction for the 
EWPOs. 

4.2 Sensitivities of “direct” tests at future colliders 

4.2.1 Future searches for sterile neutrinos produced in Z boson decays 

The estimated number of Z bosons produced by the future lepton colliders are 10® at 
the ILC [H] and 10^^ at CEPC [l5]. For the FCC-ee we use 10^^ produced Z bosons as 

‘‘We thank M. Ruan at this point for support with CEPC machine parameters. 
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Figure 7: Estimated sensitivities on the active-sterile mixing parameters at 90% confidence level from the 
EWPOs, using the uncertainties given in tab. [fi] For the CEPC, we use the fourth column labeled CEPC*. 
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Figure 8: Estimated sensitivities at 95% confidence level of a search for heavy neutral leptons at the Z pole, 
analogous to the one performed at LEP-I by Delphi [53] . For model parameters in the lighter grey area, the 
heavy neutrinos may decay more than 10 meters away from the primary vertex and thus become invisible 
to the detector (cf. |47p. 


discussed in [l6]. In analogy to the Delphi analysis discussed in sec. 13.2.11 we estimate the 
bound on achievable with a larger sample of Z bosons by scaling eq. (|53p . if no deviation 
from the SM should be found: 


(l-/i|)^(2 + ^|) 


10 "^^ 

FCC-ee 

10-9 

CEPC 

10-7 

ILC 


(62) 


where [iz = Mjmz- We show the resulting sensitivity for the future lepton colliders, based 
on this estimate, in fig. |8l The region in parameter space, where the decays of the sterile 
neutrinos take place more than ten meters away from the primary vertex, are shown by the 
grey area in the figure. We used the formula from ref. [U]. We note that even stronger 
bounds could be possible for M close to the grey shaded region, from searches for displaced 
vertices, as discussed for the FCC-ee in [46] . 


4.2.2 Searches for sterile neutrinos in M final states at 250 GeV 

One of the central aspects of a future lepton collider is the precise measurement of the 
Higgs boson properties. The dominant Higgs production mechanism is Higgs-strahlung, at 
center of mass energies around 250 GeV, which implies that large quantities of W pairs will 
be produced as a byproduct. Analogous to the analysis of the LEP-H data in sec. 13.2.21 a 
future measurement of the cross section of the process e"*“e“ ^ can be used to constrain 
sterile neutrino properties. 

The estimated W boson yield of the ILC, CEPC and ECC-ee is shown in tab. [T] where 
we include Aleph for comparison. At present, no official estimates for the systematic un¬ 
certainties are available, but discussions in the working groups are ongoing. We therefore 
only consider the statistical uncertainty for calculating our estimate. We stress that the 
estimated constraints from W boson measurements should therefore be taken with cau¬ 
tion. They rather correspond to a maximally reachable sensitivity. The estimates for the 
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Aleph 

ILC 

CEPC 

FCC-ee 

#IT’s prod. 

10 ^ 

10 ^ 

10 » 

2 X 10^ 

<5stat. on 

10 "^ 

3x10-^ 

10 -^ 

7 X 10"^ 


Table 7: Expected number of W bosons produced at the considered future colliders and at Aleph for 
comparison. We assume that the production takes place at = 250 GeV, apart from Aleph, where = 
161 to 209 GeV. Jstat. denotes the statistical uncertainty on the measurement of o^ww^ 4 e^ which we use for 
our estimates. 



— Aleph 
ILC 

-- CEPC 

— FCC-ee 



Figure 9: Estimated sensitivity on the sterile neutrino parameters from e'’'e —>■ 41 at the WW threshold 

and beyond. The black line shows the Aleph constraint for comparison. 


constraints will be updated as soon as official forecasts for the systematic uncertainties are 
available. Given this warning, we present the estimated sensitivities on the active-sterile 
neutrino mixing parameter |0e| (and lnuj) from eq. (|55p in fig. [H 

4.2.3 Sensitivities of future measurements of Higgs boson branching ratios 

The sensitivity to the Higgs boson properties at 240-250 GeV can be significantly im¬ 
proved at a future lepton collider, compared to the LHC. The currently available estimated 
precision of the ILC, CEPC and FCC-ee are shown in the tab. [HI 

For estimating the future sensitivity to sterile neutrino properties, we consider the Higgs 
branching with the highest precision, namely in tab.jHJ The results are shown in 

fig. dni A comparison of tab. [8] with hg. which includes the currently allowed branching 
ratios into heavy and light neutrinos, shows that with the increased precision the future 
measurements are indeed sensitive to deviations caused by sterile neutrinos, despite the 
already strong constraints from precision data and from the Delphi experiment. 

Furthermore, the branching ratio Hr/i^invisible also provides a promising channel for 
sterile neutrino searches, if the sensitivity for the branching ratio could reach 0.1%. In this 
case it would be complementary and comparable in sensitivity to 


21 

















Branching ratio 

ILC 

CEPC 

FCC-ee 

Bvh^ww 

6.4 

1.3 

0.9 

Brh^zz 

19 

5.1 

3.1 

BVh^'y'y 

35 

8 

3.0 

^'^e+e-^h+$T 

11 . 0 * 

3.8 

2.2 


Table 8: Estimated precision for the measurement of the Higgs boson branching ratios at future lepton 
colliders, for one year of running. The numbers are in percent, and taken from refs. [SHIS]. 

Estimated value obtained from the FCC-ee estimate rescaled with the ILC luminosity. 




Figure 10: Estimated sensitivities on the sterile neutrino properties from the decays of the Higgs boson to 
W bosons, which is the Higgs decay channel most sensitive to heavy neutrinos at future lepton colliders, 
assuming 10 years of data taking. The black line denotes the present bounds from the LHC coming from 
h —>■ 77 . 

4.2.4 Sensitivities of future measurements of e+e“ 

Another sensitive channel to search for sterile neutrino signals is the process electron- 
positron to Higgs boson plus missing transverse energy^ The SM background is given by 
e“’'e“ —> Z* —)• Z h with the subsequent decay Z —>• uv. The Higgs boson is tagged via two 
6 -jets. In the presence of sterile neutrinos, a light and a heavy neutrino can be produced 
by Z and W exchange as discussed above. The sterile neutrino then decays into a light 
neutrino and a Higgs boson, thus contributing to This process is relevant 

for M > ruH- The estimated future experimental precision for measuring this branching 
ratio is included in tab. [ 8 ] and the estimated sensitivities for the sterile neutrino properties 
are displayed in hg. [TTl 

®We would like to thank W. Murray making us aware of this possibility. 
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Figure 11: Estimated sensitivities on the sterile neutrino properties from e^e —^ ft + 


4.3 Expected future sensitivity to sterile neutrino properties from charged 
lepton flavour violation and neutrino oscillation experiments 

The current bound on the charged lepton flavour violating decay —)• ey from the MEG 
experiment [18] plays an important role in the present “indirect” probes of sterile neutrino 
properties, which we presented as results of a global fit in section 13.11 The MEG bound 
constrains the product \0e0fj\ < 1.5 x — 0.3F{xm), with F{xm) defined in eq. dHj). 

In the global fit, with a non-zero best fit value for l^el; h drives the tight constraint on |0^|. 

Euture tests of the /i —>• 3e branching ratio from Mu3e [H] and MUSIG |39l[50| and 
tests of atomic conversion rate of /r —)• e from Mu2e m and GOMET |52] have estimated 
sensitivities of order 10“^®, which can considerably improve the constraints on \9e9^\. The 
sensitivity of PRISM/PRIME [53| and a Mu2e upgrade [54] may even reach 2 x 10“^®. We 
can translate this into a sensitivity up to \9e9^\ < 3.6 x 10 “’^/y^1 — 0.3F{xm)- 

Eurthermore, the sensitivity to the branching ratio for the lepton flavour violating rare 
tau decay r —>■ ey is expected to improve to 10“® at SuperKEKB [55], which would improve 
the sensitivity to the product |0e^r| to |0e^r| < 1-5 x 10“^. The role of searches for charged 
lepton flavour violation in future global fits will depend crucially on whether deviations 
from the SM are found or not, for instance on whether the best fit value for \9(,\ will remain 
non-zero. We note that charged lepton flavour violation can also be tested at the Z pole 
at future lepton colliders, which has recently been studied for the case of the FGG-ee in 
ref. [56] . 

Finally, we would like to remark that the processes we considered in this study are 
only sensitive to the moduli \9a\ of the active-sterile mixing angles (or equivalently to the 
moduli of the Yukawa couplings \yua\)- Sensitivity to the phases of the parameters could be 
achieved in neutrino oscillation experiments, as discussed in the effective theory framework 
MUV in [571158]. 
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Direct searches 

^ Delphi (Z pole search) @2cr: IYo P . ®^=2iil®aP 

^ LHC (Higgs decays') @lcr: |y„ p ,0^=2^|e„p 

Opal (e+e“ ^ 4 leptons) @1(T: |y|=|yel, 0^=|0ep 



Other (global fit) 

— Iyl=lyvj,©"=ieep 
—■ lyl= lyvj, 0-=iej- 

— lyl=lyvj,0"=ierp 


Figure 12: Summary of present constraints on sterile neutrino properties. The LHC constraint comes from 
h — >■ 77 . 


5 Discussion and Conclusions 

We summarize the present constraints and possible future sensitivities on sterile neutrino 
parameters in the minimal symmetry protected type-I seesaw scenario (cf. section 2) in 
figs. [12] and [TSl The scenario allows for a natural realization of sterile neutrino masses 
M around the EW scale with large, even 0(1) Yukawa couplings. It is minimal but also 
sufficiently general so that a distinction between the Yukawa couplings is possible. The 
sterile neutrino parameters of the scenario, to be tested experimentally, are thus M and 

, or equivalently M and the active-sterile mixing angles Oa (with the mappings between 
the parameters summarised in table [1]). We like to argue that it is an interesting benchmark 
scenario for evaluating the present and also the future experimental sensitivities to sterile 
neutrino properties. 

Regarding the present constraints summarized in fig. [T^ we note that due to the mod¬ 
ification of the Fermi constant measured from // decays at low energies, the EWPOs also 
provide strong constraints for smaller M. In agreement with [2], we find a non-zero best fit 
value for \9e\ at 90% Bayesian confidence level. Nevertheless, to be conservative, we rather 
present our results only as constraints here. Fig. [12] also shows that the measurements of 
the Higgs branching ratios at the LHC are sensitive to decays into sterile neutrinos (only) 
in a small range around mz- 

The estimated future sensitivities are shown in fig. [13] for the ILC, CEPC and for the 
FCC-ee (TLEP). The sensitivities are qualitatively similar, however considering the current 
proposals the FCC-ee ist the most sensitive. It is interesting to note that a strong Higgs 
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program automatically leads to increased sensitivity for heavy neutrino searches around the 
IVW threshold as a byproduct. The search methods shown in the figure are sensitive to 
different (combinations of) sterile neutrino parameters. 

Provided that heavy neutrinos with mass in the considered range and sufficiently large 
Yukawa couplings exist, one would expect to obtain signals in various processes, which could 
then be used to discriminate between the active-sterile mixings |0 q| and measure/constrain 
the heavy neutrinos’ mass. Furthermore, in addition to the processes considered here, other 
searches such as for instance the search for displaced vertices as recently studied in [36] , can 
provide additional complementary information. We emphasize that the future sensitivities 
for direct searches presented here are first estimates only and a careful evaluation of the 
expected future systematic uncertainties is required for more robust forecasts. 

We also like to note that the possible sensitivity of direct searches at the Z pole at the 
FCC-ee, and, to a lesser extent also at the CEPC, are closing in on the Yukawa couplings 
from type-I seesaw models “without protective symmetry” where the expectation for the 
Yukawa couplings follows from the relation niiy ~ with nii, < 0.2 eV. The Yukawa 

couplings and the active-sterile mixing parameters, respectively, are included as dashed 
black lines in the panels of fig. [13] for comparison. Finally, for M much larger than the EW 
scale, only the indirect constraints remain and the sensitivities become independent of M 
and can alternatively be studied in the effective theory framework as done recently in [2]. 
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■» Precision constraints: lyNlyv^L 0^=l^rp 
■■■■ "Unprotected” type-I seesaw 


Figure 13: Summary of estimated sensitivities to sterile neutrino properties at future colliders. 
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